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ABSTRACT 

We present new Spitzer observations of the phase variation of the hot Jupiter HD 189733b in the 
MIPS 24 /im bandpass, spanning the same part of the planet's orbit as our previous observations in 
the IRAC 8 /im bandpass (Knutson et al. 2007). We find that the minimum hemisphere-averaged flux 
from the planet in this bandpass is 76 ± 3% of the maximum flux; this corresponds to minimum and 
maximum hemisphere-averaged brightness temperatures of 984 ± 48 K and 1220 ± 47 K, respectively. 
The planet reaches its maximum flux at an orbital phase of 0.396±0.022, corresponding to a hot region 
shifted 20 — 30° east of the substellar point. Because tidally locked hot Jupiters would have enormous 
day- night temperature differences in the absence of winds, the small amplitude of the observed phase 
variation indicates that the planet's atmosphere efflciently transports thermal energy from the day 
side to the night side at the 24 /im photosphere, leading to modest day-night temperature differ- 
ences. The similarities between the 8 and 24 fj,m phase curves for HD 189733b lead us to conclude 
that the circulation on this planet behaves in a fundamentally similar fashion across the range of 
pressures sensed by these two wavelengths. One-dimensional radiative transfer models indicate that 
the 8 fj,m band should probe pressures 2 — 3 times greater than at 24 /im, although the uncertain 
methane abundance complicates the interpretation. If these two bandpasses do probe different pres- 
sures, it would indicate that the temperature varies only weakly between the two sensed depths, and 
hence that the atmosphere is not convective at these altitudes. We also present an analysis of the 
possible contribution of star spots to the time series at both 8 and 24 /xm based on near-simultaneous 
ground-based observations and additional Spitzer observations. Accounting for the effects of these 
spots results in a slightly warmer night-side temperature for the planet in both bandpasses, but does 
not otherwise affect our conclusions. 

Subject headings: binaries: eclipsing — infrared: stars — planetary systems — stars: individual 
(HD189733) — techniques: photometric 



1. INTRODUCTION 

We currently know of more than 30 transiting plane- 
tary systems, of which the majority are gas-giant plan- 
ets orbiting extremely close (<0.05 A.U.) to their par- 
ent stars^. These planets, known as "hot Jupiters", re- 
ceive > 10,000 times more radiation from their stars than 
Jupiter does from the Sun, heating them to tempera- 
tures as high as 2000 K (Harrington et al. 2007). Most 
of these planets are expected to be tidally locked, with 
permanent day and night sides. As a result of this intense 
and highly asymmetric irradiation and their presumably 
slower rotation rates, the atmospheric dynamics of these 
planets are expected to differ significantly from those of 
the gas giant planets in the solar system. 

One of the fundamental questions regarding these plan- 
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ets is what fraction, if any, of the energy absorbed by 
the perpetually-illuminated day side is transferred to the 
night side. The answer depends on the relative sizes of 
the radiative and advective time scales and may vary 
from planet to planet depending on the specific prop- 
erties of the atmosphere. Circulation models for these 
planets (for a recent review see Showman et al. 2007) 
predict a range of possibilities, with day- night tempera- 
ture differences ranging as high as 500 — 1000 K (Show- 
man & Guillot 2002; Cho et al. 2003, 2008; Burkert et al. 
2005; Cooper & Showman 2005, 2006; Langton & Laugh- 
hn 2007; Dobbs-Dixon & Lin 2008; Showman et al. 2008). 
The form of this circulation also varies, with some models 
predicting a quasi-steady state pattern consisting of one 
or several equatorial bands of winds circling the planet 
(Showman & Guillot 2002; Burkert et al. 2005; Cooper & 
Showman 2005, 2006; Langton & Laughhn 2007; Show- 
man et al. 2008; Dobbs-Dixon & Lin 2008), and others 
predicting the formation of more complicated structures 
such as polar vortices whose positions may vary over time 
(Cho et al. 2003, 2008). 

By observing the changes in the planet's thermal emis- 
sion as a function of orbital phase, we can directly de- 
termine the day-night temperature difference for these 
hot Jupiters. Harrington et al. (2006) reported the first 
detection of these phase variations for the non-transiting 
planet v Andromedae b at 24 /.tm. If one makes reason- 
able assumptions about the predicted size and tempera- 
ture of V Andromedae b based on its mass and distance 
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from its star, the large size of the observed phase vari- 
ation imphes a large day-night temperature difference 
and correspondingly inefficient thermal homogenization 
between the day and night sides. Cowan et al. (2007) 
made similar 8 observations of three other systems, 
HD 209458, HD 179949, and 5f Peg, of which only HD 
209458 is a transiting system. They report a detection for 
the non-transiting system HD 179949, implying a large 
day-night temperature difference similar to that of v An- 
dromedae b. 

Both the observations by Harrington et al. (2006) and 
Cowan et al. (2007) are sparsely sampled, consisting of 
a series of brief visits with the Spitzer Space Telescope 
spread out over the planet's orbit, and thus do not place 
any strong constraints on the timing of minima or max- 
ima in the light curve. The best-constrained determi- 
nation of a phase curve for a hot Jupiter comes from 
Knutson et al. (2007), hereafter known as Paper I, where 
we present 8 jim observations of the transiting planet 
HD 189733b. This planet has a mass of 1.14±0.06 Mj„p 
(Bouchy et al. 2005; Torres et al. 2008) and a radius of 
1.138±0.027 Rjup (Bakos et al. 2006b; Winn et al. 2007a; 
Knutson et al. 2007; Pont et al. 2007, 2008; Torres et al. 
2008), and orbits a KOV primary with a F-band mag- 
nitude of 7.67 (Bouchy et al. 2005). At infrared wave- 
lengths it is the brightest known star with a transiting 
planet, and its favorable planet /star radius ratio makes 
it ideal for a variety of detailed measurements (see, for 
example, Deming et al. 2006; Grillmair et al. 2007; Knut- 
son et al. 2007; Tinetti et al. 2007b; Pont et al. 2008; 
Charbonneau et al. 2008; Swain et al. 2008). On UT 
2006 Oct. 28/29 we observed this planet continuously in 
the Spitzer IRAC 8 fim bandpass for 33 hours, spanning 
slightly more than half of its orbit. The high cadence 
of these data made it possible for us to fit the resulting 
phase curve with an "orange-slice" model for the planet 
consisting of twelve longitudinal strips of constant bright- 
ness. The small size of the observed phase variation ar- 
gued for highly efficient thermal homogenization between 
the two hemispheres, in contrast to the large day-night 
temperature differences inferred for v Andromedae b and 
HD 179949b. 

Although this longitudinal temperature map provides 
a wealth of information about the circulation within 
HD 189733b's atmosphere, its interpretation is compli- 
cated by the fact that the altitude of the atmospheric 
layer corresponding to the derived map depends on the 
atmospheric opacity at 8 /im. For wavelengths where 
the opacity is low the effective photosphere of the planet 
is located deep in the atmosphere, where the pressures 
and temperatures are correspondingly higher. In their 
dynamical models Cooper & Showman (2005), Dobbs- 
Dixon & Lin (2008), and Showman et al. (2008) show 
that temperatures likely become increasingly homoge- 
nized at these higher pressures, as the radiative time 
scale increases relative to the advective time scale. Thus 
we would expect that observations of the same planet 
at different wavelengths might show varying brightness 
contrasts between the day and night sides, depending on 
how deep into the atmosphere we are looking at each 
wavelength. This also means that there is an inherent 
difficulty in comparing observations of one planet {v An- 
dromedae b) at 24 /zm with observations of other planets 
(HD 179949b and HD 189733b) at 8 /xm. 



In this paper we present new observations of the phase 
variation of HD 189733b at 24 /im, spanning the same 
part of the planet's orbit as the previous observations 
at 8 /zm described in Paper I. These observations al- 
low us to compare directly the properties of this planet's 
atmosphere at different wavelengths, and to search for 
wavelength-dependent differences that might indicate the 
relative opacities and corresponding depth of the planet's 
photosphere at each wavelength. Models by Cooper & 
Showman (2005), for example, predict that the light 
curves for hot Jupiters should vary in specific ways with 
increasing depth, with the day-night brightness contrast 
decreasing and the hot region on the day side advected 
increasingly far to the east. Because these data were 
taken a year after our 8 /um observations, it is possible 
that dynamic weather patterns may have substantially 
altered the shape of the planet's phase curve (Cho et al. 
2003, 2008; Rauscher et al. 2008). If the observed light 
curve at 24 /im shows substantially different features, it 
would provide evidence for changing circulation patterns 
in the planet's atmosphere. These observations also al- 
low for a direct comparison between v Andromedae b and 
HD 189733b, providing a definitive answer to the ques- 
tion of whether the inferred day-night temperature dif- 
ferences for these two planets indicate fundamental dif- 
ferences in their atmospheres or are simply the result of 
the differing opacities in these two bandpasses. 

2. OBSERVATIONS AND ANALYSIS 

We obtained 10,104 images of HD 189733 on UT 2007 
October 25/26 using the Spitzer MIPS 24 /im array 
(Werner et al. 2004; Rieke et al. 2004) with a 10 s in- 
tegration time. Our observations spanned 35.5 hours, 
beginning 4.3 hours before the start of the transit and 
ending 2.8 hours after the end of the secondary eclipse. 
There were two interruptions for data downloads, occur- 
ring approximately 1/3 and 2/3 of the way through the 
observations. 

The standard MIPS observing sequence dithers the 
target through 14 positions on the detector array (see 
§8.2.1.2.1 of the Spitzer Observer's Manual), cychng 
several times through seven vertically-offset scan-mirror 
(chop) positions on the left side of the array with a final 
observation at the starting position and then a nod that 
places the star on the right-hand side of the array before 
repeating the same sequence. There are small differences 
in the apparent sensitivity at each position, and as a re- 
sult we elect to treat each position as an independent 
data set in our analysis. We discard a single position 
entirely, corresponding to the uppermost left position on 
the array. This position falls within a few pixels of the 
position of a bright star, HD 350998, observed during 
the second nod position when HD 189733 is located on 
the right-hand side of the array. HD 350998 is almost 
ten times brighter than our target in this bandpass, and 
as a result it is saturated in these images. This satu- 
ration may produce undesirable effects that carry over 
into other pointing positions, and indeed we find that 
the time scries for HD 189733b corresponding to this po- 
sition differs noticeably from the other positions. 

In order to calculate the flux in each image, we first 
estimate the sky background from a 41 x 41 pixel box 
centered on the position of the star, excluding the pixels 
in a central 13 x 13 pixel region that includes both the 
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Fig. 1. — Phase variation observed for HD 189733b by Spitzer 
in the MIPS 24 bandpass, with transit and secondary echpse 
visible. The data are binned in 6.3 minute intervals. In this Figure 
and Figure 2, the stellar flux as measured at the center of the 
secondary eclipse has been normalized to unity (long-dashed line). 
Panels a, b and c show the same data, but in b and c the y axis is 
expanded to show the scale of the observed variation. The bin size 
is also increased singificantly in c to allow for a better comparison 
with the model fit. The solid line in b and c is the phase curve for 
the best-fit four slice model in Figure 5, and the short-dashed line 
in b shows the expected change in the star's flux as a result of the 
rotational modulation in the visibility of star spots over the period 
of the observations. The fiux from the planet after accounting for 
these spots would be the difference between the solid and short- 
dashed lines. 



star and its M dwarf companion(Bakos et al. 2006a). We 
iteratively trim all pixels more than 3(t away from the 
median value for this subarray, make a histogram of the 
remaining pixel values, and fit a Gaussian function to the 
central region of this histogram. 

We note that the measured background drops by 1.5% 
immediately after the telescope nod moving the star to 
a new set of positions (the pattern steps through posi- 
tions 1 — 7 for 5 — 6 cycles, then nods and does the same 
for positions 8 — 14). Deming et al. (2005) noticed the 
same feature in their 24 ^m observations of HD 209458. 
In their data the measured flux from the star decreased 
along with the background flux, and they chose to take 
the ratio of the stellar to the background flux in order 
to remove this periodic drop. We find no evidence of 
a corresponding decrease in the measured flux from HD 
189733 in the images with 1.5% lower background fluxes, 
therefore there is no need to apply the correction used 
by Deming et al. (2005). In either case, if we exclude 
these images from our final binned time series we obtain 
indistinguishable results. 

We also see no evidence for the presence of the "detec- 
tor ramp" effect (Charbonneau et al. 2005; Deming et al. 
2006; Knutson et al. 2007; Harrington et al. 2007; Knut- 
son et al. 2008; Charbonneau et al. 2008) in these data. 
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Fig. 2. — Phase variation observed for HD 189733b by Spitzer 
in the IRAC 8 fira bandpass, with transit and secondary eclipse 
visible. The data are binned in 3.8 minute intervals. These data 
were originally published in Knutson et al. (2007). The stellar 
fiux as measured at the center of the secondary eclipse has been 
normalized to unity (long-dashed line). See Figure 1 for a full 
description of the plotted quantities. 



This detector effect, which has been well-characterized in 
the IRAC 8 ^m channel, causes the effective gain (and 
thus the measured flux) in individual pixels to increase 
over time. The shape of the ramp depends on the il- 
lumination level of the individual pixel, with the most 
strongly illuminated pixels (> 250 MJy Sr~^) converg- 
ing asymptotically to a constant value within the first 
two hours of observations, and the measured flux in the 
lowest-illumination pixels increasing linearly over time. 
This effect is particularly problematic for observations of 
phase variations, as the presence of this ramp can mimic 
a real rise in flux due to the planet's phase curve. In 
Paper I we found that the correction for this effect in- 
creased the uncertainty in our estimate of HD 189733b's 
8 /im night-side flux, which is based on data near the 
beginning of the observations, by a factor of five. The 
IRAC 8 ^m and MIPS 24 ^m arrays are both Si: As de- 
tectors, and thus might be expected to behave similarly. 
The IRS 16 /im array, which is also a Si:As detector, cer- 
tainly does have a detector ramp (Deming et al. 2006). It 
is not clear why the 24 fun data lacks this ramp, but we 
can immediately eliminate the shifting position of the 
star on the array as the explanation. Although IRAC 
observations of eclipses typically use a single pointing 
(Charbonneau et al. 2005; Knutson et al. 2007, 2008; 
Charbonneau et al. 2008), Deming et al. (2006) nodded 
the position of the star in their IRS 16 ^m observations 
of HD 189733, which still shows a strong ramp. Fur- 
thermore, both IRAC 8 fim and IRS 16 /j,m images also 
exhibit an even larger detector ramp in the measured 
values for the background fluxes over the period of the 
observations, independent of whether those images are 
nodded or not (Deming et al. 2006; Knutson et al. 2007, 
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Fig. 3.- This plot shows the transit a of HD 189733b observed 
by Spitzer in the MIPS 24 fim bandpass with best-fit curve over- 
plotted including timing offsets. Residuals from this fit are plotted 
in b. Data are binned in 6.3 minute intervals with error bars deter- 
mined by the RMS vaj-iation in each bin divided by the square root 
of the number of points in the bin, and the out-of-transit points 
are normalized to one. 



Fig. 4. This plot shows the secondary eclipse a of HD 189733b 
observed by Spitzer in the MIPS 24 fim bandpass with best-fit 
curve overplotted including timing offsets. Residuals from this fit 
are plotted in b. Data are binned in 6.3 minute intervals with 
error bars determined by the RMS variation in each bin divided 
by the square root of the number of points in the bin, and the 
out-of-eclipse points are normalized to one. 



2008; Charbonneau et al. 2008). We sec no evidence for a 
ramp in the sky background in our 24 /xm MIPS images, 
indicating that this effect is in fact absent from these 
images. This is consistent with the conclusions of Dem- 
ing et al. (2005), who found no evidence for a ramp in 
their 6-hour MIPS observations. It has been suggested 
that the higher background flux in the MIPS array may 
be saturating out this effect, but a comparison of MIPS 
24 /xm and IRS 16 /im peak up images indicates that 
the background flux in the 24 /im array in electrons per 
second is approximately twice that in the 16 /im array, 
which is not enough to explain the distinct behaviors of 
these two arrays. The MIPS array is run at a higher bias 
voltage than IRS; this voltage is related to the detective 
quantum efficiency of the array, but as we do not fully 
understand the origin of the ramp it is difficult to say 
whether this might explain the differing behaviors of the 
MIPS and IRS photometry. 

We estimate the flux from the star in each MIPS image 
as follows: first we subtract the background, then we fit 
the remaining flux in a circular region centered on the 
position of the star with a model point spread function. 
We use a circular region (rounded to the nearest integer 
pixel) with a radius of six pixels and a MIPS model point 
spread function for a 5000 K point source^" for these flts. 
This six-pixel radius is large enough to encompass the 
flrst Airy ring, and we find that increasing or decreasing 
the radius by one pixel does not affect the final time 
series. We fix the position of this region for all of our fits 
at a given pointing position (meaning the boundaries of 
the region used for our fits do not shift with the position 
of the star, which varies by less than half a pixel at each 
pointing position over the period of our observations). 

To fit the observed point spread function, we interpo- 
late our model to 100 times the resolution of the MIPS 
array and then rebin with the psf centered at the desired 

Available at http://ssc.spitzer.caltech.edu/mips/psf.html 



position, which we allow to vary in our fits. This allows 
us to fit for the x and y position of the star to a resolution 
of 1/lOOth of a pixel. The scatter on the final fitted posi- 
tions is typically ±0.05 pixels, five times larger than the 
resolution in our fits, so this is a reasonable choice. We 
also fit for a constant scaling factor corresponding to the 
total flux, and use the error arrays generated by the stan- 
dard Spitzer pipeline to determine the relative weighting 
for individual pixels. We note that the M-dwarf com- 
panion to HD 189733 is included within our subarray; 
we give zero statistical weight to the values within a 3 
X 3 pixel box centered on the position of the companion 
in our fits. This companion is located at a distance of 
11" from HD 189733 (Bakos et al. 2006a), which places 
it on the outer edge of the first Airy ring for our target 
star. Its fiux is only l/30th that of HD 189733 in this 
bandpass, thus a 3-pixel box is more than sufficient to 
eliminate any contribution from the companion. 

We ffag bad pixels marked by the Spitzer pipeline 
in our subarray and give them zero weight in our fits. 
To find transient hot pixels, we collect the entire set of 
702 — 840 subarray images at a given pointing position, 
and calculate the median value and standard deviation 
at each individual pixel position. We then step through 
the subarrays and mark outliers more than 3(7 away from 
the median value for that pixel position as bad pixels in 
that image. We find that 81% of our images have one or 
fewer bad pixels, and 99% have less than five bad pix- 
els in the aperture used for our fits, which contains 113 
pixels in total. This process reduces the number of large 
outliers in the final time series, although it docs not elim- 
inate such outliers completely. We found that increasing 
our threshold for bad pixels to 4c7 and then lOa outliers 
produced comparable results with an increasing number 
of large outliers in the final time series. 

After producing a time series for each pointing posi- 
tion, we iteratively select and trim outliers greater than 
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TABLE 1 
Best-Fit Eclipse Depths and Times 



Eclipse Depth Rpianet/^Star Center of Transit (HJD) O— C (s)'' 

8.0 /im Transit^ 2.387 ± 0.006% 0.1545 ±0.0002 2454037.61196 ± 0.00007 -9 ± 6 (±14) 

24 /im Transit 2.396 ± 0.027% 0.1548 ± 0.0009 2454399.24000 ± 0.00019 4 ±16 (±11) 

8.0 iim Secondary Eclipse* 0.338 ± 0.006% 2454038.72294 ± 0.00027 116 ± 23 (±6)°'<^ 

24 /(in SccondarA' EcIii>RO 0.536 ± 0.027% 21.54100.3.50.33 ± 0.00093 65 ± 80 (±11)= 



'^Knutson ot al. (2007) 

^This column gives the observed transit time minus the transit time calculated using the cphemeris derived in §2.1 
from a fit to 27 published transits of HD 189733b and the 24 /^m transit from this paper. The uncertainties arc set to 
the uncertainty in the observed transit time, while the values in parenthesis give the uncertainty in the predicted time. 
The total uncertainty in the O— C values is the sum of these two values. 

*^Predicted secondary eclipse times are defined as Tc ± 0.5P + 30 s, where the additional 30 s delay accounts for the 
light travel time in the HD 189733 system (Loeb 2005). 

*^The predicted time for this eclipse is calculated using the determined from the 8 jLim transit instead of the value 
quoted in §2.1, as this produces a more accurate prediction. 



3(7 to remove any remaining points affected by transient 
hot pixels. Wc include data spanning the transit and sec- 
ondary eclipse in this iterative trimming process, but we 
first divide the time series by the best-fit transit and sec- 
ondary eclipse light curves as determined in §2.1 before 
selecting 3(t outliers. After trimming these outliers from 
the original time series (including eclipses) we divide the 
trimmed time series by its median value and combine all 
of the pointing positions into a single time series consist- 
ing of 9,243 points, 91% of the original total. See Figure 
1 for the resulting 24 fim light curve and Figure 2 for 
the comparable 8 fxm light curve from Paper I. We set 
the uncertainties for each individual point equal to the 
standard deviation of this combined time series after the 
end of the secondary eclipse. 

2.1. Fitting the Eclipses 

We fit the transit and secondary eclipse using a Markov 
Chain Monte Carlo (MCMC) method (see, for example 
Ford 2005; Winn et al. 2007b) with 10^ steps. We ini- 
tialize the chain with the best-fit parameters determined 
from a minimization routine and add small random 
perturbations to these values to ensure that the chain 
explores the correct region of parameter space. We cal- 
culate our transit and secondary eclipse light curves us- 
ing the equations from Mandel & Agol (2002) for the case 
with no limb-darkening. Beaulieu et al. (2008) found that 
accounting for the effects of limb- darkening increased the 
resulting best-fit transit depth for HD 189733b by 0.02% 
in the 5.8 fiia IRAC bandpasses; we would expect that 
the effects of limb-darkening would be much smaller in 
the MIPS 24 /xm bandpass. As a test, we repeat our 
transit fit with a single linear limb- darkening coefficient 
as an additional free parameter, and find that the best-fit 
transit depth increases by 0.017% or 0.6 a. Our best-fit 
value for the limb-darkening coefficient is 0.022 ± 0.077, 
indicating that our data are entirely consistent with a 
limb-darkening coefficient of zero. As a result, we elect 
to fix the limb-darkening to zero in all of our subsequent 
fits. 

Our free parameters in the fit include a constant scal- 
ing factor, the transit time, and the transit depth. Af- 
ter running the chain, we search for the point where the 
value first falls below the median of all the values in 
the chain (i.e. where the code had first found an excellent 
fit), and discard all steps up to that point. We take the 



median of the remaining distribution as our best-fit pa- 
rameter, with errors calculated as the symmetric range 
about the median containing 68% of the points in the 
distribution. The distribution of values was very close 
to symmetric in all cases, and there did not appear to 
be any strong correlations between variables. Figures 3 
and 4 show the binned data with best-fit transit and sec- 
ondary eclipse curves overplotted. Best-fit eclipse depths 
and times are given in Table 1. As a check, we also re- 
peated the same fits using a standard downhill simplex 
X^ minimization routine, and obtained equivalent best-fit 
parameters. 

Because the 24 /im transit is not as well-constrained 
as the 8 /um observations described in Paper I, we set 
the inclination to its best-fit value from the fit to the 
8 ^m transit. Wc then fit for the depth of the transit at 
24 /xm, which is proportional to the square of the ratio 
of the planetary and stellar radii, and the transit time. 
We find that this 24 /zm transit depth differs from the 

8 /Ltm transit depth from Paper I by 0.3<7. We perform 
the same fit for the secondary eclipse, allowing both the 
depth and timing to vary independently and set other 
parameters to their best-fit values from Paper I. We find 
a relative depth of 0.536±0.027%. This is consistent with 
the previous value of 0.598±0.038% from a 24 nm eclipse 
observed in 2005 (Charbonneau et al. 2008), at a level of 
1.3cr. 

There have been several new high-precision observa- 
tions of transits of HD 189733b published in the past 
year, including HST ACS observations (Pont et al. 2007) 
and the Spitzer 8 /im transit from Paper I. Rather 
than using the previously published cphemeris from 
Winn et al. (2007a), which docs not include these re- 
cent observations, we derive a new cphemeris from a 
fit to all 27 previously published transits (Bakos et al. 
2006b; Winn et al. 2007a; Knutson et al. 2007; Pont 
et al. 2007). This new cphemeris has a central tran- 
sit time Tc = 2454399.23990 ± 0.00017 HJD and a pe- 
riod P = 2.21857578 ± 0.00000080 days. Using this 
new cphemeris we find that the 24 fim transit occurs 

9 ± 22 s later than predicted. Next we repeat this fit 
including the 24 /xm transit, and find a central transit 
time Tc = 2454399.23995 ± 0.00013 HJD and a period 
P = 2.21857597 ± 0.00000060 days (±52 ms). Using 
this cphemeris, wc find that the secondary eclipse occurs 
64 ± 81 s later than the predicted time, which is defined 
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TABLE 2 

Comparison of the Minimum and Maximum 

Plaxet-Star Fu'x Ratios 



Parameter 


8 ixm 


24 nm 


F - 


0.219 ± 0.024% 


0.416 ± 0.027% 


^max 


0.342 ± 0.006% 


0.550 ± 0.0027% 




64 ± 7% 


76 ± 3% 


^ rmn,corr 


0.261 ± 0.025% 


0.443 ± 0.027% 


^mirijCorr / ^max 


76 ± 7% 


81 ± 3% 



^This gives the minimum planet/star flux ratio after sub- 
tracting the estimated contributions from the star spots. 
The maximum planet /star flux ratios are not affected by 
these spots. 



as Tc + 0.5P + 30 s (the additional 30 s delay accounts 

for the light travel time between the planet and star, as 
calculated from Loeb 2005). Our timing precision is not 
sufficient to confirm or reject the f 20 ± 24 s delay in the 
time of the secondary eclipse reported in Paper I. 

2.2. Fitting the Phase Curve 

We fit the observed phase variation with an "orange- 
slice" model of the planet consisting of iVsUccs longitudi- 
nal slices with a uniform intensity in both longitude and 
latitude (Cowan fc Agol 2008). We find that smoothing 
this step function does not significantly change the result- 
ing light curve, provided the total flux from each slice and 
its brightness-weighted longitude are unchanged. The 
slices arc centered at (pojcjio + 4>o + 2A^, . . ., where 
< (/^o < a free parameter in the fit. The model 

therefore has A^siicos + 1 free parameters. The phase offset 
is necessary because it is possible to determine diQ to bet- 
ter than ±A(/)/2 since the projected area of a slice peaks 
when the center of the slice is facing the Earth. 

We test a series of models with either 2, 3, 4, or 6 
slices. The best-fit model parameters and the associated 
la uncertainties are determined using a MCMC method 
as described in §2.1. We initialize this chain using the 
best-fit values from the Levenberg-Marquardt fit. The 
reduced is fairly insensitive to changes in the number 
of slices, but the uncertainty in the intensity of each in- 
dividual slice increases with A^giiccs- We elect to use the 
four-slice model fit for our final analysis, as this repre- 
sents the best compromise between the degree of spatial 
resolution and the uncertainties in the flux from each 
individual slice. 

For all of these models, we find that the brightest re- 
gion on the planet is located to the east of the sub- 
stellar point, consistent with a peak in the integrated 
phase function occurring before the start of the sec- 
ondary eclipse. For our two-slice model the brightest 
slice is centered 31 ± 9° to the east of the substellar 
point, and in the four-slice model it is centered 22 ±8° to 
the east. The peak in the integrated light curve for 
the four-slice model, corresponding to the maximum 
hemisphere- averaged flux from the planet, occurs at an 
orbital phase of 0.396 ± 0.022 or a central meridian lon- 
gitude 37 ± 8° east of the substellar point (note that 
the longitude of the brightest slice and the longitude of 
the brightest hemisphere are not the same for asymmet- 
ric longitudinal flux distributions). With a difference of 
2.2a, this is mildly inconsistent with the 8 fiva. phase 
curve, which reaches its maximum value at an orbital 



TABLE 3 

Comparison of the Minimum and Maximum 
Hemisphere-Averaged Brightness 
Temperatures 



Parameter 


8 //m'' 


24 ^m 


^Tnax 


1258 ± 11 K 


1220 ± 47 K 


T ■ 

^ mm 


1011 ±51 K 


984 ± 48 K 


T —T- 
^ -max ^ mzn 


247 ± 51 K 


236 ± 48 K 


rr> b 
rmn,corr 


1098 ± 51 K 


1032 ± 48 K 


T" T" ■ b 

-^Tnax ^rmn,corr 


160 ± 51 K 


188 ± 48 K 



^These temperatures are higher than the values 
published in Paper I because we use a Kurucz at- 
mosphere model specific to HD 189733 (available 
at http://www.kurucz.cfa.harvard.edu/stars) to de- 
termine the flux from the star rather than interpo- 
lating from a grid of atmosphere models. 

^This gives the minimum hemisphere-averaged 
brightness temperatures after subtracting the esti- 
mated contributions from the star spots (the maxi- 
mum temperature estimates are not affected by these 
spots). 

phase of 0.456 ± 0.017 or a central meridian longitude 
16 ± 6° east of the substellar point. 

Because the region of the planet located 90° east of the 
antistellar point was visible only briefly at the beginning 
and end of the observations, it is not as well constrained 
by these data as the other regions of the planet. For the 
24 jim data, our best-fit models indicate that the cold- 
est region on the planet may be located anywhere from 
30° east (from a two-slice fit) to 90° east (from a four- 
slice fit) of the antistellar point. This would correspond 
to a minimum in the planet's phase curve that occurred 
before the start of our observations, which explains why 
its location is poorly constrained by these fits. There is 
a local minimum in the phase curve for the best-fit four- 
slice model occurring at an orbital phase of 0.1 or a cen- 
tral meridian longitude 40° west of the antistellar point, 
but this feature is not statistically signiflcant and we are 
unable to confirm or exclude the existence of a minimum 
comparable to the one observed in the 8 /zm data. 

Based on the four-slice fit, we conclude that the ob- 
served phase curve has a maximum in fiux of 1.00550 ± 
0.00010 where the stellar fiux as measured at the center 
of the secondary eclipse has been normalized to unity. 
The minimum flux in this fit is 1.00349 ± 0.00038, but 
this involves an extrapolation to a time prior to the 
start of our observations. We use a simpler two-slice 
model to estimate the; niinimuni flux and find a value of 
1.00416 ± 0.00011, which is consistent with the median 
flux value prior to the start of the transit. When we 
estimate the brightness temperatures we must include 
the uncertainty in the depth of the secondary eclipse, 
as this measurenicint dciterniines the baseline flux contri- 
bution from the star. This increases the uncertainties 
in both the minimum and maximum flux estimates to 
±0.0027. Taking the difference of the maximum four- 
slice and minimum two-slice fluxes, we find an increase 
of 0.133 ± 0.015% in the measiired flux from this system. 
This corresponds to a night-side fiux from the planet that 
is 76 ± 3% of the day-side flux, where the uncertainty in- 
cludes the propagated uncertainties in the minimum and 
maximum fluxes as well as the depth of the secondary 
eclipse, which is used to determine the planet-star flux 
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ratio. This is consistent at a level of 1.6<7 with the re- 
sults of our previous 8 /xm observations, which found 
that the night-side flux was 64 ± 7% of the day-side flux. 
The uncertainties are comparable for the 24 /xm light 
curve despite the increased scatter in the data because 
the correction for the detector ramp in the 8 /im band- 
pass adds significant uncertainty to the 8 ^m night-side 
flux (see §2 and Paper I for more detailed explanations 
of the 8 /im ramp correction) . The 24 /im night-side flux 
matches the prediction by Barman (2008), who used the 
planet's day-side broadband emission spectrum (Char- 
bonneau et al. 2008) and 8 /im phase curve to predict a 
night/day flux ratio of 75% at 24 fim. 

In order to facilitate comparisons between the model 
fits at 8 and 24 fim, we repeat our four-slice fit using the 
8 /im phase curve from Paper I. The resulting four-slice 
fit is plotted together with the 24 /im model in Figure 5, 
and the corresponding integrated phase curve is shown 
in Figure 2. We also calculate the minimum and maxi- 
mum hemisphere-averaged brightness temperatures (see 
Table 3) corresponding to the minimum and maximum 
fluxes in the 8 and 24 /im light curves from Table 2. 
To determine these brightness temperatures, we use a 
Kurucz atmosphere model for the star^^ (Kurucz 1979, 
1994, 2005) and assume a Planck function for the planet. 
We take the ratio of these two functions and integrate 
over the IRAC 8 /im and MIPS 24 /im bandpasses to 
determine the planet-star flux ratio in each bandpass, 
and solve for the temperature that matches the flux ra- 
tio observed in each bandpass. Note that our estimate 
of the planet's brightness temperature is dependent on 
our choice of the atmosphere model for the star; in Pa- 
per I we interpolated from a grid of models, whereas in 
this paper we use a model specific to HD 189733. This 
new model has a slightly higher effective temperature and 
correspondingly higher 8 /im flux, and as a result our es- 
timates of the brightness temperatures for the planet in 
the 8 /im bandpass are higher than the values given in 
Paper I. 

2.3. Effects of Starspots 

HD 189733 is an active KO star (Bouchy et al. 2005), 
which has been observed to vary by ±1.5% at visible 
wavelengths (Winn et al. 2007a; Henry & Winn 2008; 
CroU et al. 2007). Because the spots on this star have an 
effective temperature approximately 1000 K cooler than 
that of the stellar photosphere (Pont et al. 2008), the am- 
plitude of these variations scales approximately as the ra- 
tio of two blackbodies, which would imply that variabil- 
ity from spots should have a much smaller amplitude at 
24 /im. However, even a small variation in the star's flux 
during the period of our observations might contribute 
significantly to the observed signal. In order to char- 
acterize the behavior of these spots, we obtained simul- 
taneous ground-based observations of HD 189733 using 
the 1.2 m telescope at the Fred Lawrence Whipple Ob- 
servatory (FLWO) and one of the Tennessee State Uni- 
versity 0.8 m automated photometric telescopes (APT) 
at Fairborn Observatory. Although the FLWO observa- 
tions span only the week surrounding our MIPS obser- 
vations, the APT observations are part of a long-term 
monitoring program (Henry & Winn 2008) and extend 

Available at http://kurucz.harvard.edu/stars.html 
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Fig. 5. — Brightness estimates for four longitudinal strips on 
the surface of the planet at 8 /im (blue triangles) and 24 fim (red 
circles), where the brightness values are given as a percentage of 
the total flux from the planet. The solid overplotted lines are spline 
interpolations of these four-slice fits, and the dashed lines show the 
same model fits after correcting for the effects of star spots. The 
horizontal dotted lines at the bottom of the plot indicate the range 
of central meridian longitudes (longitudes viewed face-on) visible 
during our IRAC (upper line) and MIPS (lower line) observations. 

over more than a year, including the times of both our 
8 and 24 /im Spitzer data. These observations allow us 
to determine the behavior of the star during these two 
periods. 

2.3.1. FLWO Photometry 

We used KeplerCam on the 1.2 m telescope at the Fred 
L. Whipple Observatory on Mt. Hopkins, Arizona to ob- 
tain Sloan g and z photometry of HD 189733 on seven 
consecutive nights beginning on UT 2007 October 22. 
This instrument has a 23'.lx23'.l field of view, allow- 
ing us to obtain photometry of a number of bright com- 
parison stars simultaneously with our target (see Winn 
et al. 2007a, for more detailed information on observa- 
tions of HD 189733 using this instrument). HD 189733 
was setting early in the evening, so these observations 
were taken during a brief period at the beginning of each 
night. We used a single master flat constructed from the 
dome flats taken over all nights to correct the images. 
Images were intentionally defocused to avoid saturating 
the array during our 5 s exposures, and we discarded any 
images where the point spread function for HD 189733 
contained saturated pixels. We used aperture photome- 
try with IRAF's PHOT task and an aperture radius of 
35 pixels or 23.5" to estimate the flux from HD 189733 in 
each image. Increasing the size of this aperture produced 
a corresponding increase in the amount of noise from the 
sky background, while decreasing it led to larger system- 
atic variations. We estimated the level of sky background 
using an annulus with a radius 35 < r < 60 pixels. 

We corrected for variations in atmospheric transmis- 
sion and instrument efficiency using a set of comparison 
stars visible in the images, iteratively discarding compar- 
ison stars with light curves that appear to differ signifi- 
cantly from the average. There were ten comparison stars 
in the final iteration, with a relative total flux three times 
greater than that of HD 189733. We estimated the mean 
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Fig. 6. — Relative flux variation for HD 189733 observed with the 
TIO APT in the Str0mgrcn y band encompassing the time of our 
8 /irn IRAC observations (upper panel). Vertical lines indicate the 
start and end times for our Spitzer observations, which occurred 
during a period of increasing stellar flux. Variations are caused by 
rotational modulation in the visibility of star spots with a rotation 
period of 11.953 ± 0.009 days (Henry & Winn 2008). The lower 
panel shows the variations observed during a several-month period 
around the 8 //m Spitzer observations phased by the rotation period 
of the star. 

flux for HD 189733 each night relative to the ensemble 
average of these ten calibrators, and set the uncertain- 
ties equal to the RMS variation in these relative fluxes 
divided by the square root of the number of images. The 
resiilting flux values in the g bandpass for each night are 
plotted in Figure 6. We obtained similar measurements 
in z but most of the images from the latter part of the 
week had saturated pixels, and the limited time coverage 
of the remaining images meant that they were not useful 
for our analysis. 

2.3.2. APT Photometry 

We obtained observations of HD 189733 in Stromgren 
b and y filters over a span of several months surround- 
ing our 8 and 24 //m Spitzer observations from an on- 
going monitoring program carried out with the TIO 0.8 
m APT at Fairborn Observatory in southern Arizona 
(Henry 1999; Eaton et al. 2003). In these observations 
the telescope nodded between HD 189733 and three com- 
parison stars of comparable or greater brightness as de- 
scribed in Henry & Winn (2008). Because we are ul- 
timately interested in the behavior of the starspots at 
longer wavelengths, we elect to use the y band pho- 
tometry for our analysis (see Figure 6). We note that 
Henry & Winn (2008) published similar ol)servations of 
HD 189733 spanning the period of our 8 /xm observations 
of this system; this allows us to determine the behavior 
of the star during this period as well. 

2.3.3. Spitzer IRAC 8 iim Photometry 

The FLWO and APT photometry indicates that the 
star is increasing in flux during both of our observa- 
tions, and we can use this information to estimate the 
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Fk;. 7. — Relative flux variation for HD 189733 observed with the 
TIO APT in the Str0mgren y band (filled circles), with the FLWO 
1.2m in the Sloan g band (filled triangles), and Spitzer IRAC 
8 /im band (open squares) around the time of our 24 /im Spitzer 
MIPS observations. Vertical lines indicate the start and end times 
for the MIPS observations, which occurred during a period of in- 
creasing stellar flux. The first and fourth IRAC 8 /um points cor- 
respond to secondary eclipses, while the second and third points 
are transits. The average over these four points has been set equal 
to one in this plot. The first and second IRAC observations are 
located during a period of maximum stellar flux and the third and 
fourth IRAC observations are located during a period of minimum 
stellar fiux. The decrease in flux over the time of these observations 
provides a direct estimate of the size of the variations induced by 
spots in the 8 /um IRAC bandpass. 

amplitude of the corresponding variations in the 8 and 
24 /im bandpasses. We scale these variations to infrared 
wavelengths using unpublished Spitzer observations of 
two transits and two secondary eclipses of HD 189733b 
in the IRAC 8 ^lm bandpass (program GO 40238, PI E. 
Agol), taken during the weeks immediately before and 
after our MIPS observations. These observations include 
a secondary eclipse on UT 2007 October 20, a transit 
observed on UT 2007 October 21, a transit observed on 
UT 2007 November 14, and a secondary eclipse observed 
on UT 2007 November 15 (see Figure 7). The first tran- 
sit / secondary eclipse pair occur during the stellar maxi- 
mum immediately preceding our MIPS observations, and 
the second transit/secondary eclipse pair occur during 
the stellar minimum a little more than three weeks later. 

We estimate the fiux from the star in individual im- 
ages by summing the flux within circular apertures with 
radii ranging from 3.4 to 7 pixels centered on the position 
of the star. For apertures smaller than 3.4 pixels there 
are flux losses correlated with the position of the star on 
the array, and apertures larger than seven pixels become 
increasingly noisy. We determine this position using a 
weighted sum of the fluxes over a 7x7 pixel box centered 
on the approximate position of the star, following the 
methods described in Knutson et al. (2007, 2008) and 
Charbonneau et al. (2008). We subtract a background 
from each image determined by fitting a Gaussian func- 
tion to a histogram of pixels in the corners of the subarray 
images where flux from the star as minimal. We do not 
apply an aperture correction, as we are only interested 
in estimating relative changes in the flux from the star 
over time. 

We determine the decrease in the flux from the star 
during these two epochs by comparing data from the two 
transits and the two secondary eclipses separately. This 
is because we expect the flux from the planet to vary as it 
moves through its orbit, but this effect cancels out if we 
compare data from the same region of the phase curve 
(i.e. immediately after the transit or immediately af- 
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ter the secondary eclipse). This assumes that the planet 
docs not experience significant weather-related variabil- 
ity, such as that suggested by Cho et al. (2003, 2008) 
and Rauscher et al. (2008). The overall consistency in 
the observed features of our 8 and 24 /xm light curves, 
which are separated in time by approximately one year, 
indicates that weather-related variability in the observed 
fluxes is probably minimal. Our 24 /xm secondary eclipse 
depth is also consistent with a 24 /xm eclipse observed in 
2005 (see §2.1 and Charbonneau et al. 2008), placing an 
additional constraint on the variability. 

We find that the median flux measured after the end of 
the transit decreased by 0.15% between the two epochs 
and the median flux after the end of the secondary eclipse 
decreased by 0.20% between the two epochs (see Figure 
7), where each of these measurements was averaged over 
apertures ranging from 3.4 to 7 pixels in radius. We use 
a range of aperture sizes in order to test whether the de- 
tector ramp described in §2, which causes the effective 
gain of the detector to increase as a function of time, is 
affecting our result here. Pixels that are more strongly 
illuminated have both a shorter characteristic time scale 
and a smaller relative amplitude for the ramp (on the or- 
der of 1% for high-illumination pixels, versus 10% for the 
lowest-illumination pixels). Thus, over the 5-6 hour time 
frame of these eclipse observations, the high-illumination 
pixels at the center of the star's point spread fimction 
have already converged to a constant value, whereas the 
lower-illumination pixels near the edge of the aperture 
are contributing the majority of the observed ramp. We 
find that the median flux after the end of the transit de- 
creases by [0.17%, 0.12%, 0.18%, 0.14%] between the two 
epochs for apertures with radu of [3.4, 3.5, 5.0, 7.0] pix- 
els, respectively. The consistency of these values over a 
range of aperture sizes indicates that the detector ramp 
described in §2 is effectively removed in this ratio. We 
have no particular reason to prefer one aperture size over 
another, and there may be other effects at work. There- 
fore, we elect to average over these four apertures. We 
use the same method to estimate the decrease in the 
star's flux over the epoch bounded by the two secondary 
eclipse observations, and find values that are similarly 
consistent over a range of apertures. 

As a second test we take the ratio of the median fiuxes 
over the entire observations, trimming only the first hour 
of data where the ramp is steepest. This further tests 
whether the detector ramp is influencing our results, as 
we have now added in more data at earlier times when 
the detector ramp is larger. We again find consistent 
results over a range of apertures, with the star's flux de- 
creasing by 0.11% on average between one transit and 
the next, and by 0.17% on average between the two sec- 
ondary eclipses. Together with the two values derived 
solely from data after the ends of the eclipses, these val- 
ues represent four independent estimates of the decrease 
in stellar flux over this period. Combining these val- 
ues, we find that the flux from the star decreases by 
0.16 ± 0.02% in the IRAC 8 ^m bandpass between these 
two epochs, while it decreases by 1.3% over the same 
period in the Stromgren y bandpass. 

2.3.4. Scaling the Starspots to Spitzer Wavelengths 

As shown in Figure 6, we have APT observations in the 
Stromgren y bandpass spanning the period of both our 



MIPS 24 yum phase variation observations on UT 2007 
Oct. 25/26 and our previous IRAC 8 fiia phase varia- 
tion observations on UT 2006 Oct. 28/29. The FLWO g 
and APT y band observations are in good agreement and 
show that the flux from the star is increasing during the 
times of our 8 and 24 ^m Spitzer observations (these are 
the grey shaded regions in Figure 6). We determine the 
increase in the star's flux during the period of our IRAC 
8 /xm observations by fitting a linear function of time to 
the phased Stromgren y data plotted in the bottom panel 
of Figure 6, beginning at the flux minimum before the 
start of our IRAC observations and ending at the flux 
maximiim shortly after the end of these observations. 
From this fit we find that the star's flux increased by 
0.0196 ±0.0002% per hour during our IRAC 8 /xm obser- 
vations in the Stromgren y bandpass. Unlike the earlier 
data from 2006 spanning the IRAC 8 /im observations, 
the variations plotted in the top panel of Figure 6 dur- 
ing the period of our MIPS observations in 2007 do not 
phase well, indicating that the properties of the spots are 
varying during this epoch. We estimate the increase in 
the star's flux during the MIPS observations by fitting 
a linear function of time to the unphased Stromgren y 
data, beginning at the flux minimum before the start of 
our MIPS observations and ending at the flux maximum 
shortly after the end of these observations. We find that 
the star's flux increased by 0.011±0.001% per hour in the 
Stromgren y bandpass during our MIPS observations. 

To determine the star's contribution to the flux varia- 
tions observed in the 8 and 24 fxm bandpasses, we must 
scale the observed changes in the Stromgren y bandpass 
to reflect the decreased contrast of these starspots rel- 
ative to the star's photosphere at infrared wavelengths. 
For the IRAC observations, this scaling is simple: as dis- 
cussed in §2.3.3, a decrease of 1.3% in y corresponded to 
a decrease of 0.16 ± 0.02% at 8 /zm. Applying the same 
scaling to the observed increase in Stromgren y during 
our 8 /im Spitzer observations, we estimate the star in- 
creased in flux by 0.0024 ± 0.0003% per hour at 8 fim. 
The total increase in flux observed in this bandpass was 
0.12 ± 0.02% over 17.6 hours. This implies that the star 
contributed 0.042 ± 0.005% to the observed increase in 
flux, one-third of the total signal. Higher cadence MOST 
data obtained several months prior to these observations 
(CroU et al. 2007) indicates that the increase in flux dur- 
ing this part of the star's rotation is effectively linear; in 
either case the small size of the star's contribution makes 
it very unlikely that the minimum and maximum in the 
8 /im light curve can be explained by the effects of star 
spots. 

To estimate the contribution of the starspots to our 
MIPS 24 fim data, we first scale the star's 0.011±0.001% 
per hour increase in Stromgren y to the equivalent value 
of 0.0013 ± 0.0002% per hour in the Spitzer 8 /im band- 
pass. Next, we estimate how the effects of the spots scale 
between the 8 and 24 /xm bandpasses. The precise scaling 
depends on the temperatures of the spots relative to the 
effective temperature of the star. Previous HST ACS 
observations of HD 189733 have established that these 
spots have temperatures between 4000-4500 K (Pont et 
al. 2008), and we use this temperature range in our anal- 
ysis. We estimate the relative decrease in flux df (A) from 
these spots as the difference between spectra from a grid 
of model atmospheres (Kurucz 1979, 1994, 2005) with a 
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temperature of 5000 K and either 4500 or 4000 K: 



df{X) = 



h (A) - /.pot (A) 
A (A) 



(1) 



The effective temperature of the star is 5050 ± 50 K 
(Bouchy et al. 2005), so a 5000 K model is a reason- 
able choice for the star. We take the weighted average of 
df (A) over the IRAC 8 /im bandpass and then over the 
MIPS 24 Ilia bandpass, and find that a 0.0013 ± 0.0002% 
per hour increase in flux at 8 /im would correspond 
to an increase of 0.0011 ± 0.0002% per hour in the 
MIPS 24 fiTci bandpass. The observed increase in flux is 
0.133±0.015% over 25 hours, so this imphes that the star 
contributes 0.027 ± 0.004% or approximately one-fifth of 
this increase. As a check we recalculate this scaling using 
Planck functions and find that we obtain indistinguish- 
able results. 

In both the 8 and 24 /im bandpasses, accounting for the 
effects of star spots results in a slightly warmer minimum 
hemisphere-averaged brightness temperature, but does 
not otherwise alter our conclusions. This is because the 
maximum hemisphere-averaged brightness temperature 
is set by the depth of the secondary eclipse (which gives 
the total flux from the planet at that point relative to 
the flux from the star) relative to the maximum in the 
phase curve, and the interval between these two events 
is relatively short. The minimum hemisphere-averaged 
brightness temperature, on the other hand, is set by the 
changes in the observed fiux over a much longer period 
of time where the effects from star spots are increasingly 
important. We give the minimum planet /star flux ratios 
for both bandpasses after correcting for the effects of 
these spots in Table 2, and the corresponding minimum 
hemisphere-averaged brightness temperatures in Table 3. 

3. DISCUSSION 
3.1. Day and Night Atmospheric Structure 

The picture that emerges from our 24 ;um observa- 
tions broadly matches the situation we previously in- 
ferred from the planet's 8 /xm light curve in Paper I. In 
radiative equilibrium, tidally locked hot Jupiters should 
exhibit day-side temperatures at the photosphere that 
exceed 1300 K and night-side temperatures as low as 
200-300 K, implying a day-night temperature difference 
exceeding 1000 K (e.g. Showman et al. 2008; Barman et 
al. 2005). In contrast to this reference state, we find that 
the planet exhibits similar day- and night-side brightness 
temperatures at 24 /im, with a night side only modestly 
colder than the day side. In fact, the difference between 
the maximum and minimum hemisphere-averaged tem- 
peratures that we infer at 24 /im, 236 ± 48 K, is indistin- 
guishable from our previously estimated 8 jim. value of 
247 ± 51 K. Accounting for the effects of star spots re- 
sults in a slightly warmer minimum hemisphere-averaged 
temperature at both wavelengths, but does not otherwise 
affect our conclusions (see Table 3). 

Relative to these values, current three-dimensional cir- 
culation models over-predict the day- night flux variations 
at both 8 and 24 /xm (Fortney et al. 2006b; Showman 
et al. 2008). These same models indicate that the day- 
night temperature difference should increase with alti- 
tude (Cooper & Showman 2005; Showman et al. 2008; 
Dobbs-Dixon & Lin 2008), but this change is modest over 



the factor of 2 — 3 variation in the pressures that arc likely 
sensed by these two bandpasses. The small size of the 
observed flux variations in both the 8 and 24 jim. light 
curves indicates that the circulation efficiently trans- 
ports thermal energy from the day side to the night side 
over the range of pressures spanned by the 8 /xm and 
24 /im photospheres, leading to moderately (though not 
completely) homogenized temperatures between the day 
side and night side. 

There is a second possible explanation for the similari- 
ties between the 8 and 24 /im light curves. If both band- 
passes sense similar atmospheric pressures then we would 
expect to see the same features in both light curves, re- 
gardless of how the atmospheric circulation varied with 
pressure. We investigate this issue by computing day- 
side average and non-irradiated night-side ID models 
of HD 189733b that assume solar metallicity, negligible 
TiO/VO opacity due to condensation, and neglect cloud 
opacity (Fortney & Marley 2007; Fortney et al. 2008). 
The pressure-temperature proflles for these models are 
shown in Figure 8. While profiles are shown in the right 
panel, we will first examine the left panel, which shows 
normalized contribution functions (c/) for the thermal 
fiux from the atmosphere (e.g. Chamberlain & Hunten 
1987; Griffith et al. 1998): 



cf{P) = B{\T) 



Ac 



dlog(P) 



(2) 



The band-averaged c/s show the fractional contribu- 
tion of various pressures to the outgoing thermal radi- 
ation in the chosen infrared bands. These were com- 
puted by calculating the c/s at 2000 wavelengths across 
the planet's spectrum, and then, at every pressure, in- 
tegrating these c/s across the IR bandpasses. Impor- 
tantly, the peak in the c/ differs from the pressure one 
would estimate from simply solving for the point where 
the actual temperature equals the brightness tempera- 
ture (Chamberlain & Hunten 1987). The c/s show that 
there is considerable overlap in the contributions of the 
flux between the various Spitzer bands. While the peak 
in 24 /im emission is from pressures 2 — 3 times greater 
than at 8 /xm, the overlap is considerable. In general the 
shape of the c/s is the same on both hemispheres, but de- 
viations can be seen in the 3.6 and 8.0 /xm bands, which 
are affected by the increased methane abundance on the 
night side. Most prominently at 8.0 /im, the enhanced 
gaseous opacity in the upper atmosphere leads to a sig- 
nificant flux contribution from lower pressure regions. If 
one does not see as deeply on the night side compared 
to the day side, this complicates the interpretation of 
temperature homogenization. 

The c/ at 24 /im peaks at lower pressures than the 
8 /xm c/ as a result of the increased water absorption 
at longer wavelengths (e.g. Fortney & Marley 2007; Bur- 
rows et al. 2008). Figure 8 shows that for both the day 
and night sides, CO is favored over CII4, but the CII4 
abundance is not negligible, and absorption bands from 
this molecule can be seen in model spectra. The simple 
night-side proflle, that of an isolated object at Teff=1500 
K, yields a night-side synthetic spectrum that matches 
the night-side photometry, as we will show in Figure 9. 
However, this high T^s model must vastly over-predict 
the temperature of the deeper atmosphere, just below 
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Fig. 8. — Left panel: Normalized contribution functions for 
ID model atmospheres of HD 189733b. The day side is plotted 
with solid curves, while the night side uses dash-dot curves. Band- 
average contribution functions are shown for various Spitzer bands 
and at K. For clarity, the 8 and 24 fim curves are normalized to 
1.0, while the 3.6 and 4.5 fim bandpasses are normalized to 0.75 
and the remaining bandpasses are normalized to 0.5. Right panel: 
Atmospheric pressure-temperature profiles for the day and night 
hemispheres. For comparison, in dashed red is shown a ID planet- 
wide average profile, which is computed like the day-side profile, 
but with 1/2 the stellar flux, due to the day/night average. The 
night side profile is that of 1500 K isolated object, which has a 
much warmer interior than that of HD 189733b, and hence the at- 
mosphere at high pressure becomes unrealistically hot. The dotted 
curve shows where CH4 and CO are equal in abundance. 

the IRAC bands' c/s. The actual night-side Tcs must be 
considerably lower than 1500 K, to avoid an energy bud- 
get problem for the planet. To further illustrate the need 
for the planet to be relatively cool at higher pressures (as 
suggested by Barman 2008), we plot a ID planet-wide 
average profile in red in Figure 8. Much of the upper 
atmosphere of both the day and night profiles is warmer 
than this ID planet-wide average model — this cannot be 
true at all pressures. The drawbacks of simple day/night 
models computed in this fashion have been discussed pre- 
viously (Burrows et al. 2006, 2008) and argue for realistic 
3D simulations of energy transport and atmospheric tem- 
perature structures. In light of these uncertainties our 
computed c/s should be regarded with care, but as we 
show below, these same atmosphere models are a good fit 
for the mid-infrared photometric data, which argues for 
the veracity of the treatment of chemistry and opacity. 

The issue of the methane abundance on the day and 
night hemispheres is an important one. Swain et al. 
(2008) recently reported a detection of methane absorp- 
tion in the transmission spectrum of HD 189733b, which 
probes the region around the day-night terminator. It 
is reasonable to expect that the cooler night side might 
have a higher methane abundance than the day side, 
but non-equilibrium carbon chemistry could alter this 
balance. Cooper & Showman (2006) find mixing time 
scales faster than the chemical conversion time scale of 
CO to methane in hot Jupiter atmospheres, which leads 
to a day/night homogenization of the methane and CO 
abundances. At this point it is too early to make def- 
inite statements regarding this issue for HD 189733b. 
Observations of the planet's day-side emission spectrum 
around the 2.2 /ini methane band as well as the compu- 
tation of spectra from 3D dynamical models (Fortney et 



al. 2006b; Showman et al. 2008) would both help to con- 
strain the methane abundance on the day- and night-side 
hemispheres. 

If an optically thick cloud deck exists at altitudes above 
the photosphere expected from purely gaseous opacity, 
both 8 and 24 /xm radiation could emanate from the same 
pressure (that of the cloud top) and hence sense the same 
temperature structure between day and night. However, 
a high opaque cloud of silicates or iron, which are the 
most likely candidates, is not expected from equilibrium 
chemistry at these modest temperatures. Observations 
of the planet's optical transmission spectrum (Pont et 
al. 2008; Redfield et al. 2008) suggest that it may have 
a haze layer that is reducing the depth of the observed 
absorption features (for a more detailed discussion of this 
effect see Fortney 2005). However, the detection of water 
and methane absorption features in the planet's near-IR 
(1.4 — 2.5 /im) transmission spectrum indicates that this 
haze must be composed of relatively small particles, and 
it is extremely unlikely that it would affect the planet's 
emission spectrum at the relatively long (> 3 fim) wave- 
lengths of our observations (Tinetti et al. 2007b; Swain 
et al. 2008) . A thick cloud layer would also lead to a fea- 
tureless infrared emission spectrum, thus the detection 
of absorption features in the planet's 3.6 — 24 fim broad- 
band day-side emission spectrum (Charbonneau et al. 
2008; Barman 2008) provides a somewhat weaker con- 
straint on the presence of an opaque cloud layer above 
the upper range of the 8 and 24 /xm contribution func- 
tions. 

If the 8 and 24 fim photosphere pressures indeed dif- 
fer, then our observations constrain the rate at which 
the temperature changes with height in this region of 
the atmosphere. The difference between the maximum 
and minimum hemisphere-averaged 8 and 24 fim bright- 
ness temperatures is 38 ± 48 K for the maximum tem- 
perature and 27 ± 70 K for the minimum temperature. 
This suggests a structure where the temperature varies 
only weakly with pressure. In contrast, over a factor of 
two increase in pressure, a convective adiabat changes 
temperature by an amount ATad = 0.7RT/cp, where 
R, T, and Cp are the specific gas constant, tempera- 
ture, and specific heat at constant pressure. For HD 
189733b, where R « 3700 J kg'^ K"!, T w 1000 K, 
and Cp w 1.3 x 10^Jkg~^K~^, this expression yields 
ATad ~ 200 K. Thus, our observations suggest that the 
temperature increases with depth more weakly than an 
adiabat and hence that the atmosphere is not convec- 
tive at these altitudes. This result is consistent with the 
pressure-temperature profiles plotted in Figure 8 and the 
predictions of previous ID radiative-equilibrium calcula- 
tions (e.g. Fortney 2005; Seager et al. 2005; Burrows et 
al. 2005; Barman et al. 2005). Three-dimensional circula- 
tion models produce similar predictions, indicating that 
temperature should increase with pressure on the night 
side but remain close to isothermal or exhibit an inver- 
sion layer on the day side (Cooper & Showman 2005, 
2006; Showman et al. 2008; Dobbs-Dixon & Lin 2008). 

3.2. Spatially Resolved Atmospheric Features 

Our 24 fim flux maps (§2.2 and Figure 5) indicate 
that the highest-flux region lies eastward of the substel- 
lar point, providing further evidence for the horizontal 
and/or vertical advection of the temperature field by jet 
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streams, waves, or other processes. The eastward phase 
shift of the high-flux region inferred here, 20 — 30° of 
longitude, is robust to model assumptions (2, 3, 4, or 6 
slices) and is furthermore consistent with our previous in- 
ferences at 8 /um (Figure 5), which detected an eastward 
offset of ~ 30° of longitude. Our models also indicate 
that the lowest-flux region lies eastward of the antistel- 
lar point, although the size of this eastward shift varies 
substantially with model assumptions (ranging from 30° 
longitude for a 2-slice map to 90° longitude for a 4-slice 
map), due to the short duration of the data available 
before the transit. 

At first glance the 8 and 24 /im maps in Figure 5 ap- 
pear to give conflicting answers for the location of the 
cold region on the planet's night side, but the lack of 
data before the transit at 8 /xm means that we cannot 
rule out the presence of a second, larger cold spot cen- 
tered ~90° west of the substellar point, similar to the one 
hinted at in the 24 /im data. The error bars in the 12-slice 
model fit plotted in Fig. 3 of Paper I indicate that the 
fiux for the slice located 90° west of the substellar point is 
only 1(7 higher than the lowest-flux slice, which is located 
150° east of the substellar point. Without the benefit of 
the additional data before the transit available to us at 
24 fxm, we cannot determine if the minimum observed in 
the planet's integrated 8 phase curve is a local min- 
imum or a global minimum. This idea is supported by 
the new phase-shifted four-slice fit to these data plotted 
in Figure 5, which shows that the 8 //m data can be fit 
consistently by a model with two minima, one located to 
the west of the substellar point and one located to the 
east. 

Several circulation models predict that the regions of 
maximum and minimum flux will be shifted to the east 
of the substellar and antistellar points, respectively, by 
amounts analogous to that suggested by our flux maps 
(Showman & Guillot 2002; Cooper & Showman 2005; 
Fortney et al. 2006b; Showman et al. 2008). These mod- 
els do not reproduce the cold spot east of the substellar 
point indicated by the 8 /zm data (this places it in the 
same hemisphere as the hot region near the substellar 
point, which is also shifted to the east), but our new 
24 fim data argue for a global minimum located to the 
west of the substellar point, which would be consistent 
with these models. Other models suggest that the flow 
may contain hot and cold vortices that migrate in longi- 
tude; depending on timing, this could cause eastward or 
westward offsets of the flux minima and maxima (Cho et 
al. 2003, 2008; Rauscher et al. 2007, 2008). If we are to 
obtain strong constraints on the location of minima and 
other features occurring in the region to the west of the 
substellar point, it will require additional observations 
spanning the missing half of this planet's phase curve. 

3.3. Matching Models to Observations 

In light of the mid-infrared day-side emission spectrum 

(Charbonneau et al. 2008) and near infrared transmission 
spectrum (Swain et al. 2008) that have recently become 
available for this planet, it is worthwhile to take a global 
view of HD 189733b and attempt to reconcile current 
atmosphere models with our estimates for the day- and 
night-side fluxes as well as these other data. In Figure 9 
we show ratio spectra, computed from the profiles shown 
in Figure 8, as updated from Fortney & Marley (2007). 



For the planet's day side, we show a model that utilizes 
equilibrium chemical abundances (black), and in red a 
model that uses a non-equilibrium CH4/CO ratio of 0.014 
(CH4 mixing ratio of 7 x 10~^) as taken directly from 
Cooper & Showman (2006) and Fortney et al. (2006b). 
This particular CH4 abundance is consistent with the 
upper limit on the CH4 abundance on the planet's limb 
of 5 X 10~^ derived by Swain et al. (2008). Both day side 
models assume the atmosphere absorbs incident flux and 
redistributes this energy evenly over the day side, with 
no energy transported to the night side, as this provides 
the best fit to the relatively high mid-IR fluxes observed 
on the day side (blue points) . It is immediately apparent 
that both day-side models are an excellent (Icr) fit to the 
data, with the exception of the IRAC 3.6 /im point and 
the Barnes et al. (2007) K-band la upper limit. 

We model the night side as an object in isolation and 
vary the effective temperature in order to match the ob- 
served fluxes at 8 and 24 fim. Burrows et al. (2006) 
and Barman (2008) have both investigated similar kinds 
of models. In this case we find that a model with 
Teff = 1500 K provides the best fit, although this model 
has an implausibly hot interior for a planet. Barman 
(2008) find very similar fits using the 8 /im data from Pa- 
per I alone, and propose that the high night-side fluxes 
might be explained by the transport of energy from the 
day side via flows below the mid IR photosphere probed 
by Spitzer. As long as such flows remained below the 
level of the mid-IR photosphere, day-side fiux originat- 
ing from higher up in the atmosphere would appear to be 
consistent with a no-recirculation model, while the night 
side would have increased emission resulting from this 
added infiux of energy. 

As suggested by Barman (2008), increased circulation 
below the level of the mid-IR photosphere would also help 
to explain the depressed day-side iiT-band flux, as the ef- 
fective photosphere in this bandpass is below that probed 
by most mid-infrared wavelengths (see Figure 8). In this 
scenario the day-side pressure-temperature profile would 
be considerably colder than that predicted from a simple 
ID model for pressures greater than ~0.1 bar. However, 
we note that emission from the X-band and the IRAC 
3.6 /im band emerges from similar pressure ranges, such 
that a depression in fiux in A'-band but an enhancement 
in the IRAC 3.6 /xm band would be difficult to explain 
with this model. Perhaps abundant day-side methane 
coifld be suppressing the K-hand flux; this bandpass 
spans a strong methane absorption band located around 
2.3 /im. Fortney et al. (2006a) modeled this planet's at- 
mosphere and showed that a reasonable methane abun- 
dance would produce observable effects in the planet's 
if-band emission. The recent detection of methane ab- 
sorption in the planet's transmission spectrum (Swain et 
al. 2008) conflrms that this molecule is indeed present 
at the day-night terminator, although these observations 
do not provide a strong constraint on the methane abun- 
dances on the day and night hemispheres. 

If we are to reconcile fully the planet's day- and night- 
side emission spectra, it will require more sophisticated 
ID models than we have yet computed, or better yet, full 
3D simulations with non-gray radiative transfer, which 
we are working towards. Such models would go a long 
way towards characterizing both the planet's energy bud- 
get and the manner and depth dependence of tempera- 
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Fig. 9. — Planet-to-star flux ratios as a function of wavelength. Photometric data are in blue and error bars are Icr. The data from 
Charbonneau et al. (2008) are squares, the IRS 16 fim value from Deming et al. (2006) is a triangle, and our data from Paper I and this 
paper are circles. The line at 2.2 fim is the upper limit from Barnes et al. (2007). In light gray is the Grillmair et al. (2007) IRS spectrum 
ratio. The transmission functions for the Spitzer bandpasses are shown as dotted curves at the bottom of the figure. Model spectral ratios 
are shown in black and red. Diamonds show the values of these models averaged over the Spitzer bandpasses. The red model (with a fixed 
methane mixing ratio of 7 X 10~®, see text) fits the data well except for the 3.6 fim point, as does the black model, which uses equilibrium 
chemistry mixing ratios. Inset: On the night side, thin circles are the uncorrected flux ratios while the thick circles show the flux ratios 
after accounting for the effects of star spots. As on the day side, models with a flxed methane mixing ratio and equilibrium mixing ratios 
both flt the data well. 



ture homogenization in its atmosphere. Based on the 
simpler models and observations described above, we 
concur with Barman (2008), who suggested that 43% of 
the energy absorbed on the day side must be emitted on 
the night side. In Barman's framework the revised night- 
side flux ratios we present in this paper, which take into 
account the effects of star spots, argue for even more ef- 
ficient redistribution. Detections of the planet's flux at 
near-infrared wavelengths (such as J, H, or K) which 
sample the peak of the planet's emission spectrum and 
probe higher atmospheric pressures than those viewed by 
Spitzer, are required in order to constrain the redistribu- 
tion efficiency in a less model-dependent manner. 

4. CONCLUSIONS 

There are several clear conclusions that emerge from 
these observations. The planet's atmosphere exhibits 
only a modest variation in the day/night brightness tem- 
peratures at 8 and 24 /im when compared to radiative- 
equilibrium predictions for highly-irradiated, tidally- 
locked planets. This implies efhcient transport of thermal 
energy from the day side to the night side by atmospheric 
winds at the level of both the 8 and 24 /im photospheres. 
The planet's 8 and 24 /im phase curves both reach a peak 
before the secondary eclipse, indicating that the hottest 
region on the day side is shifted 20 — 30° east of the sub- 
stellar point at the location of both photospheres and 
providing additional evidence for the horizontal and/or 
vertical advection of the temperature field in the planet's 
atmosphere. The similarities between the phase curves 
at 8 and 24 fim suggests that either both wavelengths 
sense similar atmospheric pressures or that the circula- 



tion behaves in a fundamentally similar fashion across 
the relatively modest (factor of 2 — 3) range in pressures 
that atmosphere models indicate are sensed here. The 
uncertain abundance of methane in the planet's atmo- 
sphere complicates this picture (Swain et al. 2008), but 
based on other data we can definitively rule out the pres- 
ence of an opaque cloud layer above the range of the 8 
and 24 /tm contribution functions that, if present, might 
affect these observations (Charbonneau et al. 2008; Bar- 
man 2008). 

Our observations at 24 /im confirm that HD 189733b's 
phase curve is fundamentally different than that of v An- 
dromedae b or HD 179949b. Although these three plan- 
ets constitute a very limited sample, this result would 
seem to point towards the existence of two distinct classes 
of hot Jupiter atmospheres, characterized by either effi- 
cient or inefficient thermal homogenization between the 
day and night sides of the planet. This is particularly 
interesting in light of recent results by Knutson et al. 
(2008) and Charbonneau et al. (2008), who character- 
ized the day-side broadband emission spectra for both 
HD 189733b and HD 209458b. These observations re- 
vealed that HD 209458b has an atmospheric temperature 
inversion with water bands in emission instead of absorp- 
tion (Burrows et al. 2007), while HD 189733b's spectrum 
is well-described by a model with no temperature inver- 
sion and water absorption bands (Barman 2008). Al- 
though this is also a very limited sample, these results 
point towards a similar division of hot Jupiters into two 
distinct classes. 

It is possible that the presence of a temperature in- 
version and the degree of thermal homogenization may 
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be connected. Burrows ct al. (2007) and Fortney ct 
al. (2008), following up on prescient earlier work by 
Hubeny et al. (2003), have both suggested that HD 
209458b's temperature inversion might be caused by 
gas- phase TiO/VO, which would have condensed out of 
HD 189733b's cooler atmosphere. In this picture, tem- 
perature inversions would be correlated with the degree 
of irradiation, with a distinct division between the two 
classes of planets set by the condensation temperature of 
TiO/VO (Fortney et al. 2008; Burrows et al. 2008). Fort- 
ney et al. (2008) also point out that planets with temper- 
ature inversions absorb more of the incident flux higher 
in their atmospheres, where the radiative time scale is 
short compared to the advective time scale. They argue 
that this would naturally lead to large day-night tem- 
perature differences for these planets, while planets with 
lower levels of irradiation and no temperature inversions 
would be more homogenized. We note that both v An- 
dromedae b and HD 179949 have higher levels of incident 
flux than HD 189733b, placing them in the same class 
as HD 209458b. Unfortunately neither of these planets 
are eclipsing, making it difficidt to check directly for the 
presence of a temperature inversion. A better test would 
be to measure the day-night temperature difference for 
HD 209458b, which clearly does have a temperature in- 
version, or HD 149026b, whose high 8 ^m brightness 
temperature (Harrington et al. 2007) strongly favors the 
presence of a temperature inversion (Fortney et al. 2006a; 
Burrows et al. 2008). We have obtained such observa- 
tions of HD 209458b, and will report on the results in a 
future paper. 

Although HD 209458b has a temperature inversion 
while HD 189733b does not, they still share a number 
of basic characteristics, including a gas-dominated struc- 
ture with minimal to no solid core. It is not at all clear 
that information on the atmospheric dynamics of these 
planets will be applicable to smaller, core-dominated 
planets such as HD 149026b and GJ 436b (Butler et al. 
2004; Sato et al. 2005; Fortney et al. 2006a; Gillon et al. 
2007; Deming et al. 2007; Torres et al. 2007; Adams et 
al. 2008; Winn et al. 2008). These planets Ukely have 
atmospheres enriched in heavy elements, perhaps by a 
factor of 10 or more. In the solar system, there is a 
clear correlation between the percentage of planet mass 
that is core, and atmospheric metallicity (Lodders 2003). 
Uranus and Neptune have a C/H ratio of 30-40 times so- 
lar while Jupiter's is only 3 times solar. The higher 
surface gravities and potentially differing atmospheric 
compositions of HD 149026b and GJ 436b may signif- 
icantly alter the nature of the circulation between their 
presumably tidally-locked day and night sides. Moreover, 
GJ 436b is likely only pseudo-synchronized, as it has an 
orbital eccentricity of 0.15 (Deming et al. 2007; Demory 
et al. 2007), which further complicates this picture. Sec- 
ondary eclipse observations may shed some light on the 
properties of these planets, but they provide only a snap- 
shot of the global properties of the day-side atmosphere. 
Observations of the phase variations of these two planets 
would provide a considerably richer source of information 
on their spatially-resolved properties; such information is 
crucial if we are to understand the nature of atmospheric 
circulation for this distinctly different class of planets. 

Circulation models for all of these planets would also 



benefit from observations at additional wavelengths. Our 
ID radiative-equilibrium calculations indicate that the 
3.6 and 24 fim bandpasses should span the widest pos- 
sible range of pressures for HD 189733b, from 0.001 — 1 
bars. The same may hold true for HD 209458b, but its 
higher day-side temperatures should lead to an increased 
mid infrared opacity (Fortney et al. 2008), which may 
shift the c/s to lower pressures. The 4.5 and 5.8 /im band 
contribution functions overlap considerably with the 8 
and 24 /xm bands for HD 189733b, so it is possible that 
little new information would be obtained from observa- 
tions in these bandpasses. However, the 3.6 /im band 
(as well as the near-infrared J, H, and K bands) probes 
deeper atmospheric pressures that are closer to the peak 
in this planet's spectral energy distribution. For IRAC 
specifically, the 3.6 and 4.5 /xm bands are particularly 
sensitive to absorption by CO, methane, and water, and 
since they lie closer to the maximum in the planet's flux, 
give a more robust measure of the day/night effective 
temperatures. Although the intrapixel sensitivity evi- 
dent in both detectors (Reach et al. 2005; Charbonneau 
et al. 2005; Morales-Calderon et al. 2006; Knutson et al. 
2008) presents a challenge for observations of phase vari- 
ations in these two bandpasses, this effect is increasingly 
well-understood and it should be possible to develop a ro- 
bust correction with a modest investment of additional 
Spitzer time. We recommend that observations of other 
planets in these bandpasses span entire planetary orbits 
if possible; this would provide an additional check on the 
increasingly large effects of star spots at these shorter 
wavelengths, as well as resolving any ambiguities in the 
longitudinal temperature distributions of these planets. 

Such time-intensive observations in the IRAC 3.6 and 
4.5 /Ltm bandpasses would be well-matched to the pro- 
posed non-cryogenic Spitzer mission. Although Spitzer 
is predicted to run out of cryogen in spring 2009, obser- 
vations in the two shortest-wavelength channels should 
continue to achieve the same sensitivity after the cryo- 
gen is exhausted. Spitzer is the only continuing observa- 
tory that currently offers the means to study the infrared 
phase curves of extrasolar planets. Beyond Spitzer, we 
must await the launch of JWST, as such observations 
would be incredibly challenging from the ground. 
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